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Abstract 

We consider a Borel sum definition of all-orders perturbation theory for 
Minkowskian QCD observables such as the R e + e - ratio, and show that both 
this perturbative component and the additional non-perturbative all-orders 
Operator Product Expansion (OPE) component can remain separately well- 
defined for all values of energy y/s, with the perturbative component domi- 
nating as s —>■ oo, and with both components contributing as s — > 0. In the 
infrared s — > limit the perturbative correction to the parton model result 
for R e + e - has an all-orders perturbation theory component which smoothly 
freezes to the value 11(0) = 2/6, where b = (33 - 2N f )/6 is the first QCD 
beta-function coefficient, with Nf flavours of massless quark. For freezing 
one requires Nf < 9. The freezing behaviour is manifested by the "Contour- 
improved" or "Analytic Perturbation Theory" (APT), in which an infinite 
subset of analytical continuation terms are resummed to all-orders. We show 
that for the Euclidean Adler-D function, D(Q 2 ), the perturbative component 
remains defined into the infrared if all the renormalon singularities are taken 
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into account, but no analogue of the APT reorganisation of perturbation 
theory is possible. We perform phenomenological comparisons of suitably 
smeared low-energy data for the R e + e - ratio, with the perturbative freezing 
predictions, and find good agreement. 
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1 Introduction 



In this paper we wish to address the question of whether QCD perturbation 
theory can be used to make predictions in the low-energy infrared regime 
where one expects non-perturbative effects to dominate. Such an extension 
of the applicability of perturbation theory, beyond the ultraviolet regime 
of Asymptotic Freedom, would obviously enable one to test QCD in new 
ways. Reorganisations of fixed-order perturbation theory exhibiting a sta- 
ble infrared freezing behaviour have previously been formulated and studied, 
these include the so-called "Analytic Perturbation Theory" (APT) approach 
initiated by Shirkov and Solovtsov in Refs.fT (for a review see Ref.j2]) , and 
the Variational Perturbation Theory (VPT) approach Our discussion 
will address the more fundamental question of how all-orders QCD pertur- 
bation theory, and the non-perturbative Operator Product Expansion (OPE) 
contribution, can remain defined in the low-energy regime. We will discover 
that this is possible for Minkowskian observables, and that the APT approach 
should be asymptotic to the all-orders perturbative result which also exhibits 
the same infrared freezing behaviour found with APT. For Euclidean quan- 
tities, however, we will find that all-orders perturbation theory only exhibits 
stable infrared behaviour if one has complete information on the perturbative 
corrections to all-orders, and that this behaviour is not in general related to 
the infrared behaviour found using APT. 

We will focus our discussion on the R e + e -(s) ratio, at cm. energy ^/s. 
This is a Minkowskian quantity derived by analytical continuation from the 
Euclidean QCD vacuum polarization function. The corrections to the parton 
model result for R e + e -(s) will consist of a perturbative part, which can be 
developed as a power series in the renormalized QCD coupling a(s)=a s (s)/n, 
and a non-perturbative part which can be developed as an OPE in powers 
of A 2 /s, the first term corresponding to the lowest dimension relevant op- 
erator, the gluon condensate, being proportional to (A 2 /s) 2 . The key point 
is that the combination of the all-orders perturbation series and OPE must 
be well-defined at all values of s, since R e + e - (s) is a physical quantity. Each 
part by itself, however, exhibits pathologies. Specifically, the perturbation 
series exhibits n! growth in the perturbative coefficients, at large-orders n. 
Attempts to define the all-orders sum of the perturbation series using a Borel 
integral run into the difficulty that there are singularities on the integration 
contour termed infrared renormalons It turns out, however, that the 
resulting ambiguity in defining the Borel integral is of the same form as am- 
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biguities in the coefficient functions involved in the OPE, and so choosing a 
particular regulation of the Borel integral (such as principal value) induces 
a corresponding definition of the coefficient functions, and the sum of the 
two components is well-defined jH |3] . There is a further crucial pathology of 
the Borel integral, which we shall refer to as the "Landau divergence" . This 
means that at a critical energy s — Sl, the Borel integral diverges. It should 
be stressed that the value of sl should not be confused with the "Landau 
pole" or "Landau ghost" in the QCD coupling a(s). The "Landau ghost" 
is completely unphysical and scheme-dependent, whereas the divergence of 
the Borel integral is completely scheme- independent 0. For Minkowskian 
quantities such as R e + e - there is an oscillatory factor in the Borel transform 
in the integrand, arising from the analytical continuation from Euclidean to 
Minkowskian, which means that the Borel integral is finite at s = sl, and 
diverges for s < sl- To go to lower energies than sl we shall show that 
one needs to modify the form of the Borel integral, the modified form now 
having singularities on the integration contour corresponding to ultraviolet 
renormalons, correspondingly to go below s = Sl one needs to resum the 
OPE to all-orders and recast it as a modified expansion in powers of s/A 2 . 
One then finds that the ambiguities in regulating this modified Borel integral, 
are of the same form as ones in the modified OPE, and for s < sl the sum 
of the two components is again well-defined. In the infrared s — > limit the 
modified OPE resulting from resummation can contain a constant term inde- 
pendent of s even though such a term is not invisible in perturbation theory, 
so both the perturbative and non-perturbative components will contribute to 
the infrared freezing limit. The oscillatory factor in the Borel integral means 
that it freezes smoothly to 2/b in the infrared, where b = (33 — 2Nf)/6 is the 
first SU(3) QCD beta-function coefficient, with Nf quark flavours. 

The arguments sketched above suggest that the all-orders perturbative 
and non-perturbative components for Minkowskian quantities such as R e + e ~ (s) 
can separately remain defined at all energies, with the perturbative part being 
dominant in the ultraviolet and both components contributing in the infrared 
limit. One can then compare all-orders perturbative predictions with data, 
having suitably smeared and averaged over resonances [6] to suppress the 
non-perturbative mass threshold effects. In practice, of course, we do not 
have exact all-orders perturbative information. We know exactly the pertur- 
bative coefficients of the corrections to the parton model result for R e + e - to 
next-next-leading order (NNLO), i.e. including terms of order oP s |7j. Clearly, 
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conventional fixed-order perturbation theory for R e + e - will not exhibit the 
freezing behaviour in the infra-red to be expected for the all-orders perturba- 
tion theory. What is required is a rearrangement of fixed-order perturbation 
theory which has freezing behaviour in the infrared. As we have discussed in 
a recent paper [H] the resummation to all-orders of the convergent subset of 
analytical continuation terms ( "large- 7r 2 " terms) , arising when the pertur- 
bative corrections to the Euclidean Adler-D(-s) function at a given order 
are continued to the Minkowskian R e + e -(s), recasts the perturbation series 
as an expansion in a set of functions A n (s) which are well-defined for all 
values of s, vanishing as s — > oo in accord with Asymptotic Freedom, and 
with all but Ai(s) vanishing in the infrared limit, with Ax(s) approaching 2/6 
to provide infrared freezing behaviour to all-orders in perturbation theory. 
This "contour-improved" perturbation theory (CIPT) approach is equivalent 
to the Analytic Perturbation Theory (APT) mentioned above pOj in the case 
of R e + e -(s). We gave explicit expressions for the functions A n (s). At the 
two-loop level these can be expressed in terms of the Lambert VT-function 
[HIIIII]- To make contact with the all-orders perturbative result represented as 
a Borel integral, we note that the CIPT/APT reorganisation of perturbation 
theory corresponds to leaving the oscillatory factor in the Borel transform 
intact whilst expanding the remaining factor as a power series. Integrating 
term-by-term then yields the functions A n (s). The presence of the oscillatory 
factor in these integrals guarantees that the A n (s) are well-defined at all ener- 
gies. The CIPT/APT series should thus be asymptotic to the Borel integral 
at both ultraviolet and infrared energies. Whilst a reorganised fixed-order 
perturbation series exhibiting stable infrared freezing behaviour is possible 
for Minkowskian quantities, we shall show that it is not possible for Euclidean 
observables such as the Adler D-function, D(Q 2 ). For Euclidean observables 
the Borel integral does not contain the oscillatory factor and so is potentially 
divergent at s = sl, although as we shall show in the so-called leading-6 
approximation [3J |3] , the divergence is cancelled if all the infrared and ultra- 
violet renormalon singularities are included, and once again perturbative and 
non-perturbative components which are separately well-defined at all ener- 
gies can be obtained. This is only possible in the leading-6 approximation, 
however. 

The plan of the paper is as follows. In Section 2 we shall describe the 
CIPT/APT reorganisation of fixed-order perturbation theory for R e + e -, re- 
viewing the results of Ref . [S] . In Section 3 we consider how for Minkowskian 
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observables one can define all-orders perturbation theory, and the all-orders 
non-perturbative OPE in such a way that each component remains well- 
defined at all energies. The link between the all-orders perturbative result 
and the reorganised CIPT/APT fixed-order perturbation theory is empha- 
sised. We then briefly consider the corresponding problem for Euclidean ob- 
servables. In Section 4 we perform some phenomenological studies in which 
we compare low energy experimental data for R e + e -(s) with the CIPT/APT 
perturbative predictions. Section 5 contains a Discussion and Conclusions. 



2 Infra-red freezing of R e+e — CIPT/APT 

We begin by defining the R e + e - ratio at cm. energy yjs, 

JWW ^»f ^ ha + dr _°° 3) = 3^(1 + KM) . (1) 

Here the Qf denote the electric charges of the different flavours of quarks, 
and TZ(s) denotes the perturbative corrections to the parton model result, 
and has a perturbation series of the form, 



K(s) = a + J2 r na n+1 ■ (2) 

n>0 

Here a=a s (fi 2 )/n is the renormalized coupling, and the coefficients r\ and r 2 
have been computed in the MS scheme with renormalization scale fi 2 = s 
[Zj. We can consider the s-dependence of 1Z(s) at NNLO, 



s d ^r = -\pfr)= - + cn + P2 n 2 ) . (3) 

Here c = (153 — 19Nf)/12b is the second universal QCD beta-function coef- 
ficient, and P2 is the NNLO effective charge beta-function coefficient [llj, an 
RS-invariant combination of r 1; r 2 and beta-function coefficients. The condi- 
tion for lZ(s) to approach the infrared limit 1Z* as s — > is for the Effective 
Charge beta-function to have a non-trivial zero, p(lZ*) = 0. At NNLO the 
condition for such a zero is p 2 < 0. Putting Nj = 2 active flavours we find 
for the NNLO RS invariant, p 2 = —9.72, so that lZ(s) apparently freezes in 
the infrared to 1Z* = 0.43. The freezing behaviour was first investigated in 
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a pioneering paper by Mattingly and Stevenson ^2] i n the context of the 
Principle of Minimal Sensitivity (PMS) approach. However, it is not obvi- 
ous that we should believe this apparent NNLO freezing result In fact 
P2 is dominated by a large b 2 7t 2 term arising from analytical continuation 
(AC) of the Euclidean Adler D(—s) function to the Minkowskian R(s), with 
P2 = 9.40 — 7r 2 6 2 /12. Similarly the N 3 LO invariant p 3 will contain the large 
AC term — 5c7r 2 6 2 /12. This suggests that in order to check freezing we need 
to resum the AC terms to all- orders. 

R e + e - is directly related to the transverse part of the correlator of two 
vector currents in the Euclidean region, 

- M'Ms) = 4tt 2 * J d 4 xe^ < 0\T[ JfM (x)U0)}\0 > , (4) 

where s = —q 2 > 0. To avoid an unspecified constant it is convenient to take 
a logarithmic derivative with respect to s and define the Adler D-function, 

D(s) = s^U(s) . (5) 

This can be represented by Eq.(l) with the perturbative corrections lZ(s) 
replaced by 

V(s) = a + ^d n a n+l . (6) 

n>0 

The Minkowskian observable lZ(s) is related to T>(—s) by analytical con- 
tinuation from Euclidean to Minkowskian. One may write the dispersion 
relation, 

RW = jl r +, - d m , (7) 

Written in this form it is clear that the "Landau pole" in the coupling a(s), 
which lies on the positive real s-axis, is not a problem, and TZ(s) will be 
defined for all s. The dispersion relation can be reformulated as an integration 
around a circular contour in the complex energy-squared s-plane |14[ 115 j . 

TZ(s) = — [ d6V(se id ) . (8) 

One should note, however, that this is only equivalent to the dispersion rela- 
tion of Eq.(7) for values of s above the "Landau pole". Expanding T>(se zd ) as 
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a power series in d=a(se 10 ), and performing the 9 integration term-by-term, 
leads to a "contour-improved" perturbation series, in which at each order an 
infinite subset of analytical continuation terms present in the conventional 
perturbation series of Eq.(2) are resummed. It is this complete analytical 
continuation that builds the freezing of 7Z(s). We shall begin by considering 
the "contour-improved" series for the simplified case of a one-loop coupling. 
The one-loop coupling will be given by 



as 



(9) 



As described above one can then obtain the "contour-improved" perturbation 
series for 7Z(s), 



K(s) = Ax(s) +J2dnA n+1 (s) 



(10) 



n=l 



where the functions A n (s) are defined by, 



A n {s) 



1 

2^ 



d9a n 



1 

2^ 



d9- 



a n {s) 



11) 



[l + tb6a(s)/2] n ' 

This is an elementary integral which can be evaluated in closed-form as [S] 
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;i2) 



We then obtain the one-loop "contour-improved" series for TZ(s) 



IZ(s) = —arctan 
no 



7rba(s] 



a 2 (s) 



l + 6 2 7T 2 a 2 (s)/4)_ 



+d 2 



a (s) 



L(l + 6 2 7r 2 a 2 (s)/4)^ 



+. 



(13) 

The first arctan term is well-known, and corresponds to resumming the 
infinite subset of analytical continuation terms in the standard perturbation 
series of Eq.(2) which are independent of the d n coefficients. Subsequent 
terms corrrespond to resumming to all-orders the infinite subset of terms in 
Eq.(2) proportional to d\, d,2, ■ ■ ., etc. In each case the resummation is conver- 
gent, provided that \a(s)\ < 2/nb. In the ultraviolet s — * oo limit the A n (s) 
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vanish as required by asymptotic freedom. In the infrared s — > limit, the 
one- loop coupling a(s) has a "Landau" singularity at s = A™-. However, 
the functions A n (s) resulting from resummation, if analytically continued, 
are well-defined for all real values of s. Ai(s) smoothly approaches from 
below the asymptotic infrared value 2/6, whilst for n > 1 the A n (s) vanish. 
Thus, as claimed, TZ(s) is asymptotic to 2/b to all-orders in perturbation 
theory. We postpone the crucial question of how to define all-orders per- 
turbation theory in the infrared region until the next Section. We should 
also note that the functions A n (s) in Eq.(12) can also be obtained by simple 
manipulation of the dispersion relation in Eq.(7), which is defined for all real 
s. This avoids the possible objection that the contour integral in Eq.(8) is 
only defined for s above the "Landau pole" . 

Beyond the simple one-loop approximation the freezing is most easily 
analysed by choosing a renormalization scheme in which the beta-function 
equation has its two-loop form, 

^ = -^«V)(l + «0' 2 )). (14) 

This corresponds to a so-called 't Hooft scheme ^Ej in which the non- 
universal beta-function coefficients are all zero. Here c = (153 — 19Nf)/12b 
is the second universal beta- function coefficient. For our purposes the key 
feature of these schemes is that the coupling can be expressed analytically 
in closed-form in terms of the Lambert W function , defined implicitly by 
W(z)exp(W(z)) = z [m EH]. One has 



2\ 



1 



a(/i 

c[i -r vv-n 

i / \ —b/c 



c\i + W-Mv))] 

e VA- M5/ 

where A.jfg is defined according to the convention of ^H] , and is related to 
the standard definition 20\ by Ajjg = (2c/b)~ c ^ b AjYg. The "—1" subscript on 
W denotes the branch of the Lambert W function required for Asymptotic 
Freedom, the nomenclature being that of Ref. |18j. Assuming a choice of 
renormalization scale /i 2 = xs , where a; is a dimensionless constant, for the 
perturbation series of T>(s) in Eq.(5), one can then expand the integrand in 
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Eq.(6) for 7Z(s) in powers of a = a(xse ie ) , which can be expressed in terms 
of the Lambert W function using Eq.(15), 

' (16) 



c[l + W{A{s)e iK6 )\ 

where 

— If \/xs\ „ —b 



The functions A n (s) in the "contour-improved" series are then given, using 
Eqs(15,16), by 

i r i r° (—i) n 

A n (s) ee —J d6a n = -j d 9 K -J-[l + W 1 {A{s)e iKe )]- n 
1 r~x (—]) n 

+ 7T- / d9 K -^-[l + W-Misy™)} ■ (18) 

Here the appropriate branches of the W function are used in the two regions 
of integration. As discussed in Refs.piITU]. by making the change of variable 
w = W(A(s)e lK0 ) we can then obtain 



(-l) n rW-i(Ms)e>«n ^ 

2iKc n TT Jwi(A(s)e-iK*) w(l + w)' 



Ms) = i^r_ I —^i ■ (19) 



Noting that W 1 {A{s)e~ iK ^) = [W^A^e^ 71 )]* , we can evaluate the ele- 
mentary integral to obtain for n = 1, 

A^s) = 2 - - -±-Im[\n(W^(A(s)e iK ^)} , (20) 

where the 2/b term is the residue of the pole at w = 0. For n > 1 we obtain 



AJs) = \' Im 



In ( : W :} {A{s)elKn l ; ^ I J2 1 



l + W-i(A(s)e iK *)J ' ^k(l + W- 1 (A(s)e iK *)) k 

(21) 

Crucially the contribution from the poles at w = and w — — 1 cancel 
exactly. Equivalent expressions have been obtained in the APT approach 
10J. Provided that b/c > 0, which will be true for Nf < 9, we find the same 
behaviour as in the one- loop case with the A n (s) vanishing in the ultraviolet 
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2 4 6 8 



sx/A m 

Figure 1: The function Ai(s) of Eq.(20) versus sx/A^. We assume N f = 2 
flavours of quark. 



limit consistent with Asymptotic Freedom, and with A n (s) vanishing in the 
infrared limit for n > 1, and Ai(s) freezing to 2/6. To the extent that the 
freezing holds to all-orders in perturbation theory it should hold irrespective 
of the choice of renormalization scheme (RS), The use of the 't Hooft scheme 
simply serves to make the freezing manifest. In Figures 1-3 we plot the 
functions Ai(s),A 2 (s) and A 3 (s), respectively, as functions of (sx/Aj^). 
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A 3 (s) 




Figure 3: As Fig.l but for A 3 (s) of Eq.(21). 

Having shown how fixed-order perturbation theory can be reorganised so 
that it exhibits well-behaved freezing behaviour in the infra-red, we turn in 
the next section to a discussion of how all-orders perturbation theory and 
the all-orders non-perturbative OPE, can be defined in such a way that they 
remain well-defined at all energies. 

3 All-orders perturbation theory and OPE 

The corrections to the Adler D function, V(Q 2 ), can be split into a pertur- 
bative part, Vpt(Q 2 ), and a non-perturbative Operator Product Expansion 
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(OPE) part, V NP (Q 2 ), 



V(Q 2 )=V PT (Q 2 )+V NP (Q 2 ) ■ (22) 

The PT component is formally just the resummed perturbation series of 
Eq.(6), 

V PT {Q 2 ) = a{Q 2 ) + J2 d n* n+ \Q 2 ) ■ (23) 

n>0 

In addition one has the non-perturbative OPE contribution, 



v NP m = y ~ n ™ ^>:r n ^^ , (24) 



C n (Q 2 ^ 2 )<O n (fx 2 ) > 



n 



where the sum is over the relevant operators O n of dimension 2n. \i denotes 
the factorization scale, and C n is the coefficient function. For the Adler D 
function the lowest dimension relevant operator is the dimension four gluon 
condensate, 

< OIG^GJTIO > • (25) 

It will be convenient to scale out the dimensionful factor A 2n from the opera- 
tor expectation value, and combine it with the coefficient function to obtain 
the overall coefficient C n (Q 2 , /i 2 ). We can then write the V NP (Q 2 ) component 
in the form, 



;2* " 



^p(Q 2 ) = E C «^J • ( 26 ) 

We have suppressed the /i 2 and Q 2 dependence of the coefficient C n . The 
coefficients are themselves series expansions in a. 

C n = Ka 5 "(fi 2 )[l+0(a)] . (27) 

Here K is an undetermined non-perturbative normalisation, and S n is related 
to the anomalous dimension of the operator concerned. 

The definition of the all-orders perturbative component in Eq.(23) needs 
care. The series has zero radius of convergence in the coupling a. A direct 
way of seeing this is to consider the large-TV/ expansion of the perturbative 
coefficient d n , 

d n = d^NJ + d^NJ- 1 + . . . + <4°1 . (28) 
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The leading large- Nf coefficient, dn\ can be computed exactly to all-orders 
since it derives from a restricted set of diagrams in which a chain of n fermion 
bubbles (renormalon chain) is inserted in the initiating quark loop. Working 
in the so-called ^/-scheme , which corresponds to MS subtraction with scale 
= e ~ 5/3 Q 2 , one finds the exact large- N f result j2Tl 1221 123] . 



d [ : ] (V) = -^( n + l)(-iJ 



16 

n 



+ ZTTT E ^(l-2~ 2m )(l-2 2m -"- 2 )C 2m+1 



n\ . (29) 



The n\ growth of coefficients means that the perturbation series is at best 
an asymptotic one. To arrive at a function to which it is asymptotic one can 
use a Borel integral representation, writing 

poo 

V PT (Q 2 ) = / dze- z/a{ - Q2) B[V]{z) . (30) 
Jo 

Here B[D](z) is the Borel transform, defined by, 

B ^)] = • ( 31 ) 

n=0 

On performing the Borel integral term-by-term one reconstructs the divergent 
formal perturbation series for T>pt- If the series for the Borel transform has 
finite radius of convergence, by analytical continuation to the whole region 
of integration one can then define the Borel Sum, provided that the Borel in- 
tegral exists. On general grounds j2U 123 one expects that in renormalisable 
field theories the Borel transform will contain branch point singularities on 
the real axis in the complex z plane, at positions z = £fc=^ corresponding to 
infrared renormalons, IRk, k = 1, 2, 3, . . ., and at z = —Zk corresponding to 
so-called ultraviolet renormalons, UVk- Here b is the first beta-function coef- 
ficient, so that for QED with Nf fermion species b = —^Nf, whilst for SU(3) 
QCD with Nf active quark flavours, b = (33 — 2Nf)/6. Thus in QED there 
are ultraviolet renormalon singularities on the positive real axis, and hence 
the Borel integral will be ambiguous. In QCD with Nf < 33/2 flavours, so 
that the theory is asymptotically free, and b > 0, there are infrared renor- 
malon singularities on the positive real axis making the Borel integral again 
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ambiguous. For both field theories all-orders perturbation theory by itself is 
not sufficient. The presence of singularities on the integration contour means 
that there is an ambiguity depending on whether the contour is taken above 
or below each singularity. It is easy to check that, taking T> in the Borel in- 
tegral of Eq.(30) to be a generic QED or QCD observable with branch point 
singularities (1 — z/\zk\)~ Jk in the Borel plane, the resulting ambiguity for 
the singularity at z = \zk\ is of the form 



where K is complex. Using the one- loop form for the coupling, a(Q 2 ) = 
2/61n(<5 2 /A 2 ), one finds that in the QCD case, 



This has exactly the same structure as a term in the OPE expansion, Eq.(26), 
and one sees that the branch point exponent 7 of the IR renormalon is re- 
lated to the anomalous dimension of the operator, with 5k = 1 — 7fc- The idea 
is that the coefficient, in particular the constant K, is ambiguous in the 
OPE because of non-logarithmic UV divergences [2H1I2Z1- This ambiguity can 
be compensated by the IR renormalon ambiguity in the PT Borel integral, 
and so regulating the Borel integral, using for instance a principal value (PV) 
prescription, induces a particular definition of the coefficient functions in the 
OPE, and the PT and OPE components are then separately well-defined. 
That this scenario works in detail can be confirmed in toy models such as 
the non-linear O(N) a-model [2H1I2E]- For the QED case the ambiguity cor- 
responds to a Q 2 /A 2 effect. So that all-orders QED perturbation theory is 
only defined if there are in addition power corrections in Q 2 . Such effects are 
only important if Q 2 ~A 2 , here A corresponds to the Landau ghost in QED, 
A 2 ~ 10 560 m 2 , with m the fermion mass. Thus in QED such power cor- 
rections can have no phenomenological consequences and can be completely 
ignored. 

Our exact information about the Borel transform, B[D](z), for the QCD 
Adler D function is restricted to the large- Nf result of Eq.(29). In QCD we 
expect large-order behaviour in perturbation theory of the form d n ~Kn' y (b/2) n n\, 
involving the QCD beta-function coefficient b = (33 — 2Nf)/6. Motivated by 



AV 



PT 



Ke-^/^a 1 -^ , 



(32) 




(33) 
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the structure of renormalon singularities in QCD one can then convert the 
Nf expansion into the so-called 6-expansion [2H1 EDI EU E2] , by substituting 
N f = (33/2 - 36) to obtain, 

d n = d^b n + d^- l) b n ~ l + . . . + 4 0) . (34) 

The leading-6 term d^=d^b n is then used to approximate d n . Since dn = 
(— 3) n <i^6 n , it is known to all-orders from the large- Nf result. This approach 
is sometimes also referred to as "Naive Nonabelianization" [29j. It can be 
motivated by considering a QCD skeleton expansion [33], and corresponds 
to simply taking the first "one-chain" term in the expansion. It does not 
include the multiple exchanges of renormalon chains needed to build the full 
asymptotic behaviour of the perturbative coefficients, and there are no firm 
guarantees as to its accuracy. The leading-6 result for the Borel transform of 
the Adler-D function in the ^/-scheme can then be obtained from Eq.(29). 

B[v i %z) = fm ± 'Mil + Ajg) f.mU) ± *M 

3=1 \ ~ Zj) v 22 1 j=3 V x zj) 

so that one sees in the leading-6 limit a set of single and double pole renor- 
malon singularities at the expected positions. The residues of the UVj poles, 
A (j) and Ai(j), are given by [30J 

= 8(-lW + 6j + 2) = 4M-D J+1 (2 J + 3) 



Because of the conformal symmetry of the vector correlator [31] the IRj 
residues, B (j) and B\{j), are directly related to the UVj ones, with B (j) = 
-M-j) and Mj) = -M-j) for j > 2. B (l) = B 1 (l) = B 1 (2) = 0, and 
Bq(2) = 1. Notice the absence of an IRi renormalon singularity. This is 
consistent with the correspondence between OPE terms and IR renormalon 
ambiguities noted above, since there is no relevant operator of dimension 2 
in the OPE. The singularity nearest the origin is then the UV\ singularity 
at z — —2/6, which generates the leading asymptotic behaviour, 

€ >J2^,,(±)\ (37) 



27 V 2 , 

We shall now consider the correction, 7Z(s), to the parton model result 
for R e + e -. This may be split into a perturbative component 1Zpt{s), and 
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an OPE component TZnp(s), analogous to Eqs.(23),(24). Inserting the Borel 
representation for T*pt of Eq.(30) into the dispersion relation of Eq.(7) one 
finds the representation 

i r~ s+ie rit r°° 
K PT (s) = — - dze- z ^B[V](z). (38) 

^J-s-ie t JO 

It will be convenient to consider the all-orders perturbative result in leading- 
b approximation to start with, in which case the coupling a(t) will have its 
one-loop form, a(t) = 2/ (&ln(t/Ay)), where we assume the V^-scheme. In this 
case the t integration is trivial and one finds, 

= r dze -^^p-BlV^] M , (39) 
Jo nbz/2 

where B[D^](z) (in the ^-scheme) is given by Eq.(35). It is now possible 
to explicitly evaluate TZ^{s) in terms of generalised exponential integral 
functions Ei(n,w), defined for Rew > by 



Ei(n,w) = j dt— . (40) 
One also needs the integral 

jTto-^^^^), (41) 

Writing the 'sin' as a sum of complex exponentials and using partial frac- 
tions one can then evaluate the contribution to lZ^{s) coming from the UV 
renormalon singularities, i.e. from the terms involving Aq(j) and A\(j) in 
Eq.(35) [HI 

-JZ PT (s)\uv = ~ b [j- - yj arctan 

2 j= °° 

+ ~ h E (MJ)M^J) + (MJ) - A 1 (j>,)0 + (2, j)) ,(42) 

where (2 — tt 2 /6 is the Riemann zeta-function, and we have defined 

<f) + {p, q)=e z " /ais) {-l) q Im[Ei{p, (l/a(s)) + t7cbz q /2)] . (43) 
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To evaluate the remaining contribution involving the IR renormalon singu- 
larities we need to regulate the integral to deal with the singularities on 
the integration contour. For simplicity we could choose to take a principal 
value prescription. We need to continue the Ei(n, w) defined for Rew > 
by Eq.(40), to Rew < 0. With the standard continuation one arrives at a 
function analytic everywhere in the cut complex w;-plane, except at w = 0; 
with a branch cut running along the negative real axis. Explicitly |3o"j 



/ \n— 1 

Eiln, w) 



(n- 1)! 



71—1 j 

— lmu - 7b + > — 
L — 4 m 



m=l 



-w) 



T ? 1 ^ x r . (44) 

/ — ' (m — n + l)m 

m=0 v 1 



with 7£ = 0.5722. . . Euler's constant. The Inw contributes the branch cut 
along the negative real w-axis. To obtain the principal value of the Borel 
integral one needs to compensate for the discontinuity across the branch cut, 
and make the replacement Ei(n,w) — > Ei(n,w) + msign{Imw). This leads 
one to introduce, analogous to Eq.(43), 



0_(p, q) = e- z " /a{s \-l) q Im[Ei(p, (-l/a(s)) - mbz q /2)\ 

irRe[((l/a(s)) + inb/2) p - 1 } . (45) 



(p-1)! 

The principal value of the IR renormalon contribution is then given by J^OJ 
,11),.,, _ 2 f 14 8^ {nba(s) \ 2B a (2) 



oo 

The perturbative component is then the sum of the UV and (regulated) IR 
contributions, 

nf T {s) = nfl(s)\ uv + n^( s )\ m 

= A arct an(^^J + A^(A o ( J )0 + (l, J ) + (Ao(j)- A(j>,)0+(2,j)) 

' 3=1 

+ ^^04l,2) + ^f;(i?o(j>4l,jO + (5o(jO + 5i(j>,)0-(2,j)) (47) 



nb Tib 

3=3 
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Note that the £2 contributions cancel, and one obtains the arctan term, which 
is the leading contribution, Ai(s), in the CIPT/APT reformulation of fixed- 
order perturbation theory. The connection between the Borel representation 
and the A n (s) will be further clarified later. 

We now turn to the infrared behaviour of the regulated Borel integral. 
In the one-loop (leading-6) case the V^-scheme coupling a(s) becomes infinite 
at s = sl=A v . The e~ z l a ^ term in the Borel integrand approaches unity 
at s = sl, but the trigonometric factor sm(nbz/2)/(iibz/2) ensures that the 
integral is defined at s = sl- For s < sl, however, a(s) becomes negative, 
and the e~ z / a ^ factor diverges at z = 00, the Borel transform in the V- 
scheme does not contain any exponential ^-dependence to compensate, so 
the Borel integral is not defined. We shall refer to this pathology of the 
Borel integral at s = sl as the "Landau divergence". It is important to 
stress that the Landau divergence is to be carefully distinguished from the 
Landau pole in the coupling. The Landau pole in the coupling depends on 
the chosen renormalization scale. At one-loop choosing an MS scale [i 2 = xs, 
the coupling a(xs) has a Landau pole at s = Aj^/x, the Borel integral of 
Eq.(39) can then be written in terms of this coupling as, 

Tlf T {s) = ^N^ [se V3f z / 2 B[I )W](,) . (48) 
Jo nbz/2 

In a general scheme the Borel transform picks up the extra factor [rre 5 / 3 ] 6 ^ 2 
multiplying the V-scheme result. The Borel integrand is scheme (x) invariant. 
The extra factor has to be taken into account when identifying where the 
integral breaks down, and one of course finds the Landau divergence to be 
at the same x-independent energy, s = sl = e 5//3 A|^ = Ay. Thus the Borel 

representation of Eq.(38) for TZp^(s) only applies for s>s L . For s < Sl the 
one-loop (V-scheme) coupling a(s) becomes negative. We can rewrite the 
perturbative expansion of 1Zpt(s) as an expansion in (— a(s)), 

n PT (s) = a(s)+r 1 a 2 (s)+r 2 a 3 (s) + ... + r n a n+1 (s) + ... 

= -[(-a(s)) - n(-a(s)) 2 + r 2 (-a(s)f + ... + (-l) n r n (-a(s)) n+1 + . 

The expansion in (— a(s)) follows from the modified Borel representation 

poo 

1Z PT {s) = - dze- z ^- a ^B[R](-z) 
Jo 
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r — oo 

= / dze z/{ - a{s)) B\K](z) . (50) 
Jo 

This modified form of Borel representation will be valid when Re(a(s)) < 
0, and involves an integration contour along the negative real axis. Thus, 
it is now the ultraviolet renormalons UVk which render the Borel integral 
ambiguous. The ambiguity in taking the contour around these singularities 

~ k 

(analogous to Eq.(33)) now involves (s/A 2 ) . Of course, it is now unclear 
how these ambiguities can cancel against the corresponding OPE ambiguities. 
The key point is that since only the sum of the PT and OPE components is 
well-defined, the Landau divergence of the Borel integral at s = sl, must be 
accompanied by a corresponding breakdown in the validity of the OPE as an 
expansion in powers of (A 2 /s), at the same energy. The idea is illustrated 
by the following toy example, where the OPE is an alternating geometric 
progression, 




At any value of s, TZnp(s) is given by the equivalent functions in the middle 
line. For s > A 2 these have a valid expansion in powers of A 2 /s, the standard 
OPE, given in the top line. For s < A 2 the standard OPE breaks down, but 
there is a valid expansion in powers of s/A 2 given in the bottom line. Thus 
for s < sl the OPE should be resummed and recast in the form, 

KNp(s) = ^C n (jp\ . (52) 

It is crucial to note that this reorganised OPE can contain a Co term which 
is independent of s, as indeed is the case in the toy example of Eq.(51). 
Of course, an analogous C term in the standard OPE in Eq.(26) is clearly 
excluded since it would violate Asymptotic Freedom, and all the terms in 
the regular OPE are perturbatively invisible. As a result TZnp(s) can have 
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a non- vanishing infrared limit, and both components can contribute to the 
infrared freezing behaviour. It should be no surprise that perturbation the- 
ory by itself cannot determine the infrared behaviour of observables, but the 
existence of a well-defined perturbative component which, as we shall claim, 
can be computed at all values of the energy using a reorganised APT version 
of fixed-order perturbation theory, is a noteworthy feature. The remaining 
terms present in this modified OPE should then be in one-to-one correspon- 
dence with the UV n renormalon singularities in the Borel transform of the 
PT component, and the PT renormalon ambiguities can cancel against cor- 
responding OPE ones, and again each component separately be well-defined. 
The modified coefficients C n will have a form analogous to Eq.(27), 

C n = Ka in (ji*)[l + 0(a)] . (53) 

The anomalous dimension is that of an operator which can be identified using 
the technique of Parisi [21]. The anomalous dimension corresponding to C\ 
for the Adler D function has been computed j3B]- The ambiguity for the 
modified Borel representation of Eq.(50), taking UVk to be a branch point 
singularity (1 — zj ' Zk) lh , is 

ARpT^Ka 1 -^ (j^j . (54) 

Comparing with Eq.(53) one finds 5k = 1 — %■ The modified Borel represen- 
tation for Ti!>pT valid for s < sl will be, 

n { p ] T {s) = - / dze- z ^- a ^B[JZ (L) ](-z) . (55) 
Jo 

This may again be written explicitly in terms of Ei(n,w) functions. One 
simply needs to change a(s) — > — a(s), Zj — > —Zj, and A\{j) —> — Ai(j), 
Bi(j) — > —Bi(j) in Eq.(47). One finds that the result of Eq.(47) is invariant 
under these changes, apart from the additional terms which we added to the 
Ei(n,w) in continuing from Rew > to Rew < 0, in order to obtain the 
principal value. In fact the PV Borel integral is not continuous at s = sl- 
Continuity is obtained if rather than the principal value we use the standard 
continuation of the Ei(n,w) defined by Eq.(44). That is we redefine 

<P4p,q)=e~ z ^ s) {-l) q Im[Ei{p,{-l/a{s)) - mbz q /2)\ . (56) 
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This simply corresponds to a different regulation of singularities. We then 
sec that Eq.(47) for 7^(s) is a function of a(s) which is well-defined at all 
energies, and freezes to 2/b in the infra-red. We note that the branch of the 
arctan changes at s = sl, so that its value smoothly changes from zero at 
s = oo to 7r at s = 0. The reformulated OPE of Eq.(52) together with the 
perturbative component determines the infrared freezing behaviour, and in 
the ultraviolet the perturbative component dominates. The key point is that 
both components can be described by functions of s which are well-defined 
at all energies. The apparent Landau divergence simply reflects the fact that 
the Borel integral and OPE series, which are used to describe the PT and 
NP components, each have a limited range of validity in s. The connec- 
tion with the CIPT/APT rearrangement of fixed-order perturbation theory 
is now clear. It is obtained by keeping the sm(7rbz/2)/(7ibz/2) term in the 
Borel transform intact, and expanding the remainder in powers of z. Or- 
dinary fixed-order perturbation theory, of course, corresponds to expanding 
the whole Borel transform in powers of z. The retention of the oscillatory 
sin factor in the Borel transform ensures that the reformulated perturbation 
theory remains defined at all energies. One then finds that for s>sl, 

. / -. f°° , _ z/a(s) sin{nbz/2) z n ~ l 

A ^ = J ^ M) lbz/2 (n-l)l ' (57) 
where the one-loop A n (s) are given by Eqs.(12). Similarly for s<s L one finds 

n{) ~J ^ nbz/2 (n-1)!- { } 

Thus the CIPT/APT fixed-order result should be an asymptotic approxima- 
tion to the Borel integral at both large and small values of s. In Fig. 4 we 
compare the all-orders leading-6 result for IZpr(s) given by Eq.(47), with the 
NNLO CIPT/APT prediction, 

TZf PT (s) = A 1 (s) + d[ L) A 2 (s) + d { 2 L) A 3 (s) . (59) 

The one-loop A n (s) are given by Eqs.(12) and as in Eq.(47) the V^-scheme 
is assumed. We assume Nf = 2 quark flavours. One sees that there is good 
agreement at all values of s/Ay. 

We now turn to the full QCD result beyond the one-loop approximation, 
and as in Section 2 it will be sufficient to consider the two-loop result since 
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one can always use an 't Hooft scheme. Consider the Borel representation 
for TZpt(s) of Eq.(38). We shall assume that, as in the leading-5 approxi- 
mation, the Borel transform B[D](z) in the V^-scheme does not contain any 
exponential dependence on z, but is simply a combination of branch point 
singularities. It is then clear that the Landau divergence occurs when the fac- 
tor e~ z / a(_s ) becomes a diverging exponential, that is when Re(l/a(—s)) < 0. 
Thus the critical energy sl is given by the condition Re(l/a(—s)) = 0. At 
one-loop level one has 



2 Va 2 



v 



+ 



mb 



(60) 



and so the condition yields s = s L = Ay, as we found before. At the two-loop 
level the situation is slightly different. Integrating the two-loop beta-function 
in Eq.(14) now gives, 



a(-s) 



+ cln 



ca(—s 



1 + ca{— s) 



-In 



A 2 



mb 



(61) 



The vanishing of Re(\/a{—s)) then corresponds to the solution of the tran- 
scendental equation 



Re < cln 



ca(— s) 
1 + ca(— s) 



2 \A 2 V 



(62) 



0.4574A y . Since the Borel 



Assuming Nf = 2 flavours one finds s = s L 
integral is scheme-invariant so must the value of sl be, in particular the 
breakdown of the Borel representation would occur in any scheme, not just 
an 't Hooft one. We can perform the t-integration in Eq.(38) in closed form, 
and arrive at the two-loop Borel representation 

' e -z/a(-s+ie) 



n PT (s) 



-2 

Tib 



f 

Jo 



-s + ie) 



(Izlin ■ cc~-' Ei | L - - - — ) B[V\{z) 

(63) 

The factor in the square bracket plays the role of the e z / a ( s ^sin(7r6z/2) / (irbz/2) 
factor in the one-loop case. It provides an oscillatory factor so that at s = s L 
the Borel representation remains defined. For s < s L one must switch to a 
modified Borel representation as in Eq.(50), writing 



n PT {s) = -^-l ^ dze-^-^B[V]{-z). 



(64) 
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Which, performing the t-integration gives 



n PT {s) 



2 



f 

Jo 



dzim 



e ~z/(~a(-s+ie)) 



- ce~ zc Ei [l,-zc + 



(—a(—s + ie)) 



B[D](-z) . (65) 



Unfortunately we cannot write down a function analogous to Eq.(47) which 
gives TZpt(s) at all energies, because we do not know B[D](z) exactly. The 
two-loop situation, however, is the same as that at one-loop. The regulated 
representation of Eq.(63) applies for s>Sl , with the corresponding stan- 
dard OPE. Below s = sl one needs the modified representation of Eq.(65) 
together with the resummed OPE recast in the form of Eq.(52). The pertur- 
bative component IZpt(s) then freezes to 2/6 in the infra-red, we can see this 
if we split B[V](-z) into (1 + (B[V](-z) - 1)). The part of the integrand 
proportional to B\D]{— z) — 1 vanishes for all z from — > oo in the infra-red 
limit. The remaining term integrates to give us Ai(s) which freezes to 2/b 
as s — > 0. The non-perturbative component TZnp(s) given by the reformu- 
lated OPE together with the perturbative component determine the infrared 
freezing behaviour. There is again a direct connection with the CIPT/APT 
reformulation of fixed-order perturbation theory. Using integration by parts 
one can show that that for s>sl 



A n (s) = — - / dzim 
Kb J 



-z/a(—s+ie) 



ce zc Ei [l,zc + 



a(—s + ie) 



v n-l 



(n-1)! ' 
(66) 

where the A n (s) correspond to the two- loop results in Eqs. (20,21). Once 
again CIPT/APT corresponds to keeping the oscillatory function in the Borel 
transform intact, and expanding the remainder in powers of z. Similarly for 
s<sl one has, 



A n (s) — — I dzim 
™ Jo 



e -z/(-a(-s+ie)) 



ce 



Ei 1, -ZC + 



(—a(—s + ie)) 



n-l 



[n 



1)! 



(67) 



Thus, as in the one-loop case, the CIPT/APT reformulation of fixed-order 
perturbation theory will be asymptotic to the Borel representations at small 
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Figure 4: 6K(s) = K>p±(s) - V>Xpt{s), at the one loop level for 2 flavours of 
quark. 



and large energies. We would like, as in Fig. 4 for the one-loop case, to com- 
pare how well the fixed-order CIPT/APT perturbation theory corresponds 
with the all-orders Borel representation. We are necessarily restricted to us- 
ing the leading-6 approximation since this is the extent of the exact all-orders 
information at our disposal. One possibility is to simply use the leading-6 
result for the Borel transform, B[D^](z), in the two-loop Borel represen- 
tation of Eq.(63). The difficulty though is that with a(—s) the two-loop 
coupling, the Borel integral is now scheme-dependent, since B[D^](z) has a 
scale dependence which exactly compensates that of the one-loop coupling. 
Using a renormalization scale /i 2 = xs our result for Hpt{s) has an unphys- 
ical x- dependence. This difficulty is exacerbated if we attempt to match the 
result to the exactly known perturbative coefficients di and d,2, which we 
could do by adding an additional contribution (d\ — d^)z + (e^ — d^)(z 2 /2) 
to the Borel transform. Thus, as has been argued elsewhere, such match- 
ing of leading-6 results to exact higher-order results yields completely ad hoc 
predictions, which may be varied at will by changing the renormalisation 
scale jS3 EE]- The resolution of this difficulty follows if one accepts that 
the standard RG- improvement of fixed-order perturbation theory is incom- 
plete, in that only a subset of RG-predictable UV logarithms involving the 
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energy scale s are resummed. Performing a complete resummation of these 
logs together with the accompanying logs involving the renormalisation scale, 
yields a scale-independent result. This Complete Renormalisation Group Im- 
provement (CORGI) approach applied to T>(s) corresponds to use of a 
renormalisation scale ft 2 = e~ 2d ^ b s, where d denotes the NLO perturbative 
correction d\ in Eq.(23), in the MS scheme with /i 2 = s. In the CORGI 
scheme we have the perturbation series, 

V(t) = a (t) + X 2 aftt) + X 3 a 4 (t) + ... + X n a n + l + . . . , (68) 

where a (t) is given by Eq.(15) with z = (—l/e)(\^i/Ao) b ^°, where Ao=e d ^ b Aj^g, 
and X n are the CORGI invariants, and only X2 is known. We can then at- 
tempt to perform the leading-6 CORGI resummation, 



^corgM = ao(t) + X 2 a 3 (t) + J2xL L) a n +1 (t) + . . . , (69) 

n>2 



so that the exactly known NNLO X 2 coefficient is included, with the remain- 
ing unknown coefficients approximated by X$ , x[ L \ . . ., the leading-6 ap- 
proximations. We stress that a (t) denotes the full CORGI coupling defined 
in Eq.(15). One can define this formal sum using the Borel representation 
of T> in Eq.(30), with the result for B[D^] in Eq.(35). The integral can 
be expressed in closed form in terms of the Exponential Integral functions 
Ei(n,w) of Eq.(40), with the result 



00 

rt L \l/a{t)) = Y>{ - e^/ a ^Ez(l,z 3 /a(t))[( Zj /a(t))(A (j) - Zj A x {j)) - Zj A x {j)\ 



3=1 

+ (A (j) - Zj A x {j))} - e-«' a ®B {2)Ei(l, -z 3 /a(t)) 
00 

+ £ i - e~ Zj/a{t) Ei(h -zs/«t))[(*iWt)){B°V) + ZjMM 

J'=3 

- (B (j) + ZjB^j))} . (70) 

To define the infra-red renormalon contribution we have assumed the stan- 
dard continuation of Ei(n,w) from Rew > to Rew < 0, defined by 
Eq.(44). In a principal value was assumed, which corresponds to adding 
—i7Tsign(Im(zj/a(t)) to the Ei(l,—Zj/a(t)) term. As we found for 7£pj.(s) 
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the principal value is not continuous at s = sl, whereas the standard con- 
tinuation is. The formal resummation in Eq.(69) then corresponds to jU], 

V&rgM = V {L) (-^ + d?>(V)) + (X 2 - X { 2 L) )a 3 (t) , (71) 

once again a (t) is the full CORGI coupling, and df\v) denotes the NLO 
leading-6 correction in the ^-scheme. Inserting T> C oRGi{t) inside the disper- 
sion relation of Eq.(7) one can then define, 

n^ ORGI {s) = i r +le dt ^omM . (72) 



2ttz „ 

This can be evaluated numerically, if we have T^coRGii 8 ^) then we can obtain 

(73) 

If we set Si to be large enough we can evaluate T^corgi( s ^) using the circular 
contour in the s- plane, as in Eq.(8). Combining this circular integral with 
the integrals above and below the real negative axis we arrive at T^corgi( s ^) 
where s 2 can be as far into the infrared as we want. The all-orders CORGI 
result can be compared with the NNLO CIPT/APT CORGI result, 

n APT (s) = A 1 (s) + X 2 A 3 (s) . (74) 

Here the A n (s) are the two-loop results of Eqs. (20,21), with A(s) = (— l/e)(y / s/A D ) 
in the CORGI scheme. Analogous to Fig. 4 we plot in Fig. 5 the comparison 
of the all-orders and NNLO APT CORGI results, Nf = 2 quark flavours are 
assumed. As in the one-loop case there is extremely close agreement at all 
values of s. For the fits to low-energy R e + e - (s) data to be presented in the 
next section, therefore, we shall use the NNLO CORGI APT result. 

Before turning to phenomenological analysis in Section 4, we conclude 
this section with a brief discussion of the situation for Euclidean observables. 
We can define the Adler D function in the Euclidean region by inverting the 
integral transform corresponding to the dispersion relation of Eq.(7). That 
is we can write, 

r°° ds 

V{Q 2 ) = Q 2 -——,TZ(s). (75) 
Jo {s + Q 2 ) 



-b/c 
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One can certainly define a Euclidean version of APT by inserting the Minkowskian 
A n (s) in the right-hand side of Eq.(75), and defining 



Ai E \Q 2 ) = Q 2 
The one-loop result would be [I] 



ds 



(s + Q 2 ) 



AJs) 



2 n 



(76) 



a[ e \q' 2 



A 2 



ln(Q 2 /A 2 ) A 2 -Q 2 



(77) 



This Euclidean APT coupling freezes in the infrared to 2/6, but this be- 
haviour is induced by the second non-perturbative contribution, which can- 
cels the forbidden tachyonic Landau pole singularity present in the first per- 
turbative term. There is now no direct connection, however, between this Eu- 
clidean APT coupling and the Borel representation for Dpj-(Q 2 ) of Eq.(30). 
Since there is now no oscillatory factor present in the Borel integral it is 
potentially divergent at s = sj> We can explicitly exhibit this divergent 
behaviour working in leading-6 approximation. The Borel integral can then 
be explicitly evaluated in terms of Ei functions as we have seen in Eq.(70). 



29 



Using Eq.(44) for the Ei function one then finds a divergent behaviour as 
s — > s L proportional to lna, 

oo 

Vpt(s) -> [J2(z*Mj) + ZjBxti)) - ^£o(2)]hia + . . . , (78) 

where the ellipsis denotes terms finite as s — > sl- However, remarkably, 
the factor in the square bracket vanishes, and the result is finite at s = sl, 
provided that all the renormalon singularities are included. The contribution 
of any individual renormalon is divergent. The cancellation follows because of 
an exact relation between the residues of IR and UV renormalons (Eq.(36)), 

zjA l (j) = -zj+ a B l (j + 3). (79) 

This results in cancellations in the sum, leaving a residual term z|I?i(3) 
which then cancels with the Z2-Bo(2) term. An analogous relation Ao(j) = 
—B Q (j + 2) has been noted in [3U|. It seems that these relations are un- 
derwritten by the conformal symmetry of the vacuum polarization function 
[3*4] . but further investigation is warranted. The above finiteness at s = sl 
means that one can obtain a T>pt{Q 2 ) component well-defined in the infrared 
by changing to the modified form of Borel representation for s < sl- One 
finds that this component becomes negative before approaching the freezing 
limit Ppt(O) = 0. Similar behaviour is found for the Gross-Llewellyn Smith 
and polarised and unpolarised Bjorken structure function sum rules, whose 
complete renormalon structure is also known in leading-6 approximation [3U| . 
Phenomenological investigations are planned [HI] ■ Comparable investigations 
in the standard APT approach have been reported in [IT]. Unfortunately, 
nothing is known about the full renormalon structure beyond leading-6 ap- 
proximation. Such knowledge would be tantamount to a full solution of the 
Schwinger-Dyson equations. Correspondingly no analogue of the APT reor- 
ganisation of fixed-order perturbation theory asymptotic to T>pt is possible 
in the Euclidean case. 

We finally note that in the case of R e + e - and D it is possible to say 
something about the separate infrared freezing behaviours of the PT and 
NP components. Arguments of spontaneous chiral symmetry breaking in the 
limit of a large number of colours jSl] imply that D(0) = 0, or equivalently 
T>(0) = —1. Furthermore according to Ref.|12] H and T> should have the 
same infrared freezing limit. This argument follows directly from Eq.(8) if 
the circular contour is shrunk to zero. These exact results then suggest that 
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Figure 6: Comparison of CORGI APT and the standard NNLO CORGI 
calculations of R e + e - (s) at low energies. 



£Vp(0) = —1 to be consistent with the leading-6 result Dpt(O) = obtained 
above. For 71 one infers that TZ^p(0) = — 1 — (2/b) to be consistent with the 
n PT (0) = 2/b result. 

4 Comparison of NNLO APT with low energy 
R e + e ~ data 

In this section we wish to compare the NNLO CORGI APT perturbative 
predictions with low energy experimental data for R e + e -. The discussion 
so far has assumed massless quarks. To include quark masses we use the 
approximate result 03] 

R e+e - (s) = 3 Q 2 f T(v f ) [i + g{v f )n{s)\ , (80) 
/ 

with the sum over all active quark flavours, i.e. those with masses m/ < 
y/s/2, and where 

Vf = (l-4mj/s) 5 , 
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T(v) = v(3-v 2 )/2 



4tt 

9(v) = -j 



IT 3 + V ( 7T 3 



2v 4 V2 4tt 



(81) 



For the theoretical predictions we shall take 7Z(s) to be the NNLO CIPT/APT 
CORGI result of Eq.(74). Starting with AgL = 216MeV for N f = 5 , cor- 
responding to the world average value a s (Mz) = 0.1172 jUj, we demand 
that 7l(s) remains continuous as we cross quark mass thresholds. This 
then determines for Nt = 4, 3, 2. We take standard values for cur- 

rent quark masses for the light quarks |44] : m u = 3.0MeV, = 6.75MeV, 
m s = 117.5MeV, and also from jH] we take the values for pole masses of the 
heavy quarks m c = 1.65GeV, and mj = 4.85GeV. The approximate result 
|Hj uses pole masses in Eq.(81), so we use pole masses where we can. Using 
these values for the quark masses and a s (Mz), we plot the resulting R e + e - (s) 
in Fig. 6. The solid line corresponds to the CORGI APT result for 7Z(s) in 
Eq.(74). The dashed curve corresponds to the standard NLO fixed-order 
CORGI result, 

n CO RGi(s) = a {s) + ( X 2 - —J a 3 (s) . (82) 

The standard fixed-order result breaks down at s = A 2 D = 0.4114GeV 2 , 
where there is a Landau pole. The APT result smoothly freezes in the infra- 
red. The dashed-dot curve shows the parton model result (i.e. assuming 
K(s) = 0). 

For a recent comprehensive review of the experimental data for R e + e - (s) 
at low energies see Ref . j3H] • The experimental data we have used comes from 
a variety of sources. From the two pion threshold up to t/s = 1.43GeV 
we use references |46j . the data from these references is given as individual 
exclusive channels which must be combined to obtain the full hadronic cross 
section. In the region between 1.43GeV and 2.0GeV we use data from |4*7] . 
[4*Hj . references [4*H] . |KH] are used in the region between 2.0GeV and 5.0GeV. 
From 5.0GeV to 7.25GeV we use [55, and from 7.25GeV to 10.52GeV we 
use [32J, jSS]- These sets of data all give the inclusive total hadronic cross 
section. Above 10.52GeV we insert the NNLO CORGI APT prediction for 
R e + e -, this is represented by the continuous line in Fig. 7. 

In order to simplify the analysis of the data we did not use overlapping 
datasets, instead where one dataset overlapped another we simply took the 
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Figure 7: Data used to compare with model, statistical errors shown only. 



better, smaller error, dataset in the region of the overlap in a/s. Errors were 
dealt with by taking each data point and calculating the effect of its statis- 
tical and its systematic error. The effect of its statistical error was added in 
quadrature with the other statistical errors. The contribution from the sys- 
tematic error was added to the other systematic errors from the same dataset, 
then the contribution from the systematic errors of each dataset were added 
in quadrature with each other and the contribution from the statistical errors. 

We also need to consider the effect of narrow resonances not included in 
the data, we employ the same approach as used in [12]. We assume that the 
narrow resonances have a relativistic Breit-Wigner form 



9 

Rres(s) = —zBuB, 

cr 



M 2 r 2 



[s - M 2 ) 2 + M 2 r 2 ' 



(83) 



where a is the QED coupling, and M, T, Bu, B^ are the mass, width, lepton 
branching ratio, and hadron branching ratio respectively. We are assuming 
a narrow resonance i.e. T is small, so we approximate the resonance with a 
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delta function 

9 MTln 9 

Rres(s) = ^B u B h MTn {s _ M2)2 , + M2T2 « - 2 B u B h MTn5(s - M 2 ). 

(84) 

The compilation of data for R e + e - is shown in Fig. 7. Narrow resonances 
are indicated by the vertical lines. Unfortunately it is not possible to directly 
compare the experimental data with the theoretical predictions. This is be- 
cause there is not a direct correspondence between the quark mass thresh- 
olds in perturbation theory and the hadronic resonances. This difficulty can 
be overcome if one employs a "smearing procedure". We shall employ the 
method proposed by Poggio, Quinn and Weinberg 0, defining the smeared 
quantity 

A y~ Re+e-(t) 
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R e+e - (s; A) — — J dt - ; 2 W A2 • (85) 



o 



(t — s) + A 



R e + e -(s) itself is related to the vacuum-polarization function Tl(s) of Eq.(4) 

by, 

2iR e+e - (s) = II(s + ie) - U(s - ie) , (86) 

that is it is the discontinuity across the cut. The smeared R e + e -(s; A) can 
be written as 

2iR e+e - (s; A) = U(s + zA) - - iA) . (87) 

If A is sufficiently large one is kept away from the cut, and is insensitive 
to the infrared singularities which occur there. If both data and theory are 
smeared they can then be compared. In this way one hopes to minimise 
the contribution of the IZnp component. One needs to choose A sufficiently 
large that resonances are averaged out. For the charm region it turns out 
that A = 3GeV 2 is a good choice, whilst for lower energies A = lGeV 2 
is adequate. In Fig.8(a) we choose A = lGeV 2 . -R e + e -(s;A) obtained 
from the data is represented by the solid line. The dashed-dot line is the 
smeared NNLO CORGI APT prediction, assuming the quark mass thresh- 
olds as above with the exception of the charm quark whose mass is taken to 
be m c = 1.35GeV for reasons which we shall shortly discuss. The dashed 
line is the parton model prediction. The shaded region denotes the error in 
the data. It is clear that in the charm region the averaging is insufficient, 
although for lower energies the agreement is extremely good. In Fig. 8(b) we 
show the corresponding plot with A = 3GeV 2 . There is now good agree- 
ment between smeared theory and experiment over the whole s range, for 
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m c = 1.35GeV. Whilst we have indicated an error band associated with the 
data, we have not indicated an error band for the theory prediction. There 
are several potential sources of error to consider. The first is the choice of 
renormalisation scale. Our viewpoint would be that the use of the CORGI 
scale corresponds to a complete resummation of ultraviolet logarithms, which 
in the process results in a cancellation of /^-dependent logarithms contained 
in the coupling a(/z 2 ) and in the perturbative coefficients. As we have argued 
elsewhere attempts to estimate a theoretical error on the perturbative 
predictions by making ad hoc changes in the renormalization scale are sim- 
ply misleading, and give no information on the importance of uncalculated 
higher-order corrections. A common approach, for instance, is to use scales 
fi 2 = xs where x is varied between x = ~ and x = 2, with x = 1 providing 
a central value. We should note, however, that were we to have used such 
a procedure it would not have led to a noticeable difference in the theory 
curves, since the APT has greatly reduced scale-dependence, as has been 
noted elsewhere [S3]. A more important uncertainty is the precise value of 
the quark masses assumed, and in particular the choice of the charm quark 
mass m c . To illustrate how this effects the results we show in Fig.8(c) the 
curves obtained if we assume m c = 1.65GeV. As can be seen the theory 
curve is now inconsistent with the data in the charm region, although for 
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lower energies where the charm quark has decoupled, the agreement is again 
good. 

The uncertainty in the mass of the charm quark is exceptionally large. 
Looking at the different references used in jH] a value m c = 1.35GeV for 
the pole mass is reasonable, and agrees well with [33] which is referenced in 
|44j . Part of the problem is the relationship between the pole mass and the 
MS mass for the charm quark, where the oP s contribution is larger than the 
a 2 contribution. Obtaining the pole mass through MS mass calculations, 
which is done in jHj, is not very satisfactory. Reference joT . which also fits 
low-energy R e + e - data, gives a pole mass of m c = 1.33 — 1.4GeV, and so the 
choice of 1.35GeV is reasonable. 

It is possible to extend the smearing to spacelike values of s. We give the 
corresponding curves for .R e + e - (s; A), with m c = 1.35GeV, over the range 
-3 < s < 1 GeV 2 in Figs. 9 (a), 9(b), for A = lGeV 2 , and A = 3GeV 2 , 
respectively. The agreement between theory and experiment is extremely 
good in both cases. 

In Fig. 10 we show _R e + e -(.s;A) in the upsilon region. The choice A = 
lOGeV 2 works quite well, we show the theory predictions for different rrib 
values. A direct comparison between theory and data which does not involve 
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smearing is possible if one evaluates the area under the R e + e -(s) data, that 
is evaluates the integral, 



J(s) = / R e+e -(t)dt, 

J 4m| 

where s lies well above the low-energy resonances in the continuum. We 
show the theory and experimental I(s) over the range 5 < y/s < 9GeV in 
Fig. 11. There is extremely good agreement. Finally we can avoid smearing 
by transforming R e + e -(s) to obtain D(Q 2 ) in the Euclidean region, using the 
dispersion relation of Eq.(75) 

D(Q 2 )=Q 2 aS R e+e -(s). (89) 

U m \ [s + Q 2 ) 

In practice we cannot integrate up to infinity so we just take the sufficiently 
large upper limit of 10 6 GeV . As noted earlier above \/s = 10.52GeV the 
NNLO CORGI APT prediction is used for the data. The theory and data 
results are shown in Figs. 12,13. There is good agreement. Our results are 
comparable to the fit obtained in [56 j , and to the results in [57J- We should 
also note that very similar plots and fits to those we have presented in this 
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Section are included in Ref.jHH], which uses instead the so-called Variational 
Perturbation Theory (VPT) approach [3]. 

5 Discussion and Conclusions 

The Analytic Perturbation Theory (APT) approach advocates the "analyti- 
zation" of the terms in standard perturbation theory so that the perturbative 
expansion is recast as an expansion in a basis of functions that have desirable 
analytic properties, in particular the absence of unphysical "Landau poles" 
in Q 2 pQ. The functions in the Euclidean and Minkowski regions are interre- 
lated by the integral transforms of Eq.(7) (V — > TZ) and Eq.(75) (1Z — > V). 
In a previous paper we pointed out the Minkowskian formulation of APT for 
the quantity R e + e - was equivalent to the all-orders resummation of a con- 
vergent subset of analytical continuation terms jHj- This reorganisation of 
fixed-order perturbation theory gives apparent infrared freezing to the limit 
2/6 to all-orders in perturbation theory, and the functions A n (s) at two-loop 
level could be written in closed form in terms of the Lambert W function. 
However, one might question whether this all-orders perturbative freezing 
has any physical relevance. It is well-known that all-orders perturbation 
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Figure 13: Same as figure 12 but viewed over a smaller range. 

theory by itself is insufficient, and that it must be complemented by the 
non-perturbative Operator Product Expansion (OPE) jUE]. It is clear that 
the OPE breaks down as s — > 0, since it is an expansion in powers of A 2 /s. 
In this paper we have shown how both the PT and the OPE components 
can remain defined in the infrared limit. Writing a Borel representation for 
the PT component one finds that it is ambiguous because of the presence 
of singularities on the integration contour, termed infrared renormalons [5]. 
These ambiguities, however, are of precisely the same form as OPE terms, 
and a regulation of the singularities in the Borel integrand induces a defi- 
nition of the OPE coefficients, allowing the two components to be defined. 
We showed that the Borel integral representation inevitably breaks down at 
a critical energy Sl which we referred to as the "Landau divergence". For 
Minkowskian quantities the Borel Transform contains an oscillatory factor 
which means that the Borel integral remains defined at s = sl- For s < sl 
one needs to switch to an alternative Borel representation, which has ambi- 
guities due to ultraviolet renormalon singularities on the integration contour. 
Correspondingly the OPE should be resummed and recast in the form of an 
expansion in powers of s/A 2 . The UV renormalon ambiguities in the Borel 
integral are then of the same form as the terms in the modified OPE, and 
regulating the modified Borel integral induces a definition of the coefficients 
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in the modified OPE, allowing both components to be defined. The modified 
Borel integral freezes to 2/6 in the infrared thanks to the presence of the 
oscillatory factor, whilst the modified OPE component will also contribute 
to the infrared freezing behaviour since resummation of the standard OPE 
can result in s-independent terms which can give a nonzero freezing limit, 
as in the toy example of Eq.(51). As we noted we did not expect to be able 
to determine the infrared behaviour from perturbation theory alone, but the 
existence of a perturbative component which can be defined using a reorgan- 
ised version of fixed-order perturbation theory at all energies is important. 
In particular the perturbative component dominates in the ultraviolet and 
may possibly provide a good approximation into the low-energy region. We 
explicitly constructed the all-orders Borel representations using the all-orders 
leading-6 approximation for TZ(s) t 3 j , and a one-loop coupling. We could 
express the Borel integral in closed form in terms of exponential integral 
functions (Eq.(47)). With the standard continuation of the Ei(n,w) func- 
tions defined by Eq.(44) the result for TZ^{s) of Eq.(47) is a function of 
s which is well-defined at all energies, freezing to 2/6 in the infrared, and 
continuous at s = s^. The two-loop Borel representation was also discussed. 
The details are similar to the one-loop case, with a modified oscillatory fac- 
tor and a shifted value of sl, the modified Borel representation again freezes 
to 2/6 in the infrared. At both one- loop and two-loops the APT modifica- 
tion of fixed-order perturbation theory corresponds to keeping the oscillatory 
factor in the Borel integrand intact, and expanding the remainder. As a re- 
sult the APT results should be asymptotic to the Borel representations at 
all energies, underwriting the validity of the all-orders perturbative freezing 
behaviour. It should be noted that we have somewhat oversimplified our 
discussion of the OPE contribution. The OPE coefficients are not constant, 
as in the toy example of Eq.(51), but are functions of a, C n (a). Each co- 
efficient will involve a perturbation series in a which is divergent with n\ 
growth of coefficients, and can be defined using a Borel representation. As 
defined by analytic continuation from the OPE for T>np to that for TZnp, 
the corresponding Borel integrands will contain the same oscillatory factors, 
enabling C n (a) to remain defined at s = sl, and for s < sl one switches 
to the modified Borel representation. We should note that the difficulty of 
uniquely extending the Borel representation for Minkowskian quantities into 
the infrared has also been discussed in Ref-jSHl, but with differing conclu- 
sions to us. A more closely related discussion concerning the significance and 
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interpretation of the Landau Pole is given in Ref . |32| . The modified Borel 
representation of Eq.(50) and the promotion of UV renormalon singularities 
to the positive axis in the Borel z-plane has also been discussed in Ref. [34J. 

Whilst the Minkowskian version of APT is underwritten by a Borel rep- 
resentation valid at all energies, this is not the case for the Euclidean version. 
There is no oscillatory factor in the integrand in the Euclidean case, and the 
Borel integral will potentially diverge as one approaches sl- However, we 
showed that working in leading-6 approximation 1>pt was finite at sl thanks 
to a cancellation between the infinite set of IR and UV renormalon residues. 
For individual renormalon singularities the Borel integral is divergent. By 
switching to the modified Borel representation one can then define a Dpy 
component which in fact freezes to zero in the infrared. This is interesting 
and similar perturbative freezing is also found for structure function sum rules 
|40j . The key point, however, is that no analogue of the Minkowskian APT 
reorganisation of fixed-order perturbation theory is possible in the Euclidean 
case, and one is restricted to the leading-6 approximation in exhibiting the 
perturbative freezing. 

In the final Section we performed fits of NNLO APT results to low energy 
R e + e - data. We needed to introduce quark masses approximately, and in or- 
der to avoid ambiguities due to the precise location of quark mass thresholds, 
and to minimise the contribution of the TZnp component, we used a smearing 
procedure. Extremely good agreement between theory and data was found. 

An obvious further application would be to use the APT approach in 
the analysis of the tau decay ratio and in particular the estimation of the 
uncertainty in a s (Mz) which such measurements imply 0121]. In Ref. |DJ this 
was estimated by comparing NNLO CIPT in the CORGI approach, with an 
all-orders resummation based on the leading-6 result. However, in fact CIPT 
for the tau decay ratio is not equivalent to the APT approach and corresponds 
to an expansion in a different basis of functions. In particular the resulting 
functions are discontinuous at s = sl- We hope to study this further in a 
future publication. 
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